ChemistiyOPEN 

*** 



Open Access 



DOI: 10.1002/open.201100015 

Halogen Bonding versus Hydrogen Bonding: A Molecular 
Orbital Perspective 

Lando P. Wolters and F. Matthias Bickelhaupt* [a] 



We have carried out extensive computational analyses of the 
structure and bonding mechanism in trihalides DX-A~ and the 
analogous hydrogen-bonded complexes DH-A~ (D, X, A = F, CI, 
Br, I) using relativistic density functional theory (DFT) at zeroth- 
order regular approximation ZORA-BP86ATZ2P. One purpose 
was to obtain a set of consistent data from which reliable 
trends in structure and stability can be inferred over a large 
range of systems. The main objective was to achieve a detailed 
understanding of the nature of halogen bonds, how they re- 



semble, and also how they differ from, the better understood 
hydrogen bonds. Thus, we present an accurate physical model 
of the halogen bond based on quantitative Kohn-Sham molec- 
ular orbital (MO) theory, energy decomposition analyses (EDA) 
and Voronoi deformation density (VDD) analyses of the charge 
distribution. It appears that the halogen bond in DX-A~ arises 
not only from classical electrostatic attraction but also receives 
substantial stabilization from HOMO-LUMO interactions be- 
tween the lone pair of A~ and the a* orbital of D-X. 



Introduction 

Hydrogen bonds are, without doubt, one of the most impor- 
tant intermolecular interactions known today. Being responsi- 
ble for the unique features of water, as well as playing a key 
role in DNA structure and replication, the importance of hydro- 
gen bonding for human life can hardly be overestimated. 
Therefore, it is no surprise that this interaction has been stud- 
ied extensively. 11,21 Halogen bonds, although discovered around 
150 years ago, 131 have received considerably less attention. 
After the first experimental studies on this phenomenon by 
Hassel, 141 a review by Bent appeared in which donor-acceptor 
interactions with halogens were discussed. 151 Early theoretical 
descriptions were given by Pimentel, 161 Mulliken, 171 and Flurry.™ 
In the last decades, there has been an increased interest in hal- 
ogen bonding, which nowadays has applications in various 
fields in chemistry, 191 such as, supramolecular, 110 " 121 biochem- 
istry 113 " 151 and inorganic chemistry. 1161 Earlier studies generally 
indicate that halogen bonds can, both in terms of practical ap- 
plications and bond strength, compete with hydrogen 
bonds. 117 " 241 

In this study, we have computationally investigated a range 
of strongly halogen-bonded trihalides DX-A~ and the analo- 
gous strongly hydrogen-bonded complexes DH-A~ (D, X, A = 
F, CI, Br, I), using relativistic density functional theory (DFT). The 
purpose of our work is twofold. Firstly, we wish to provide 
a set of consistent structural and energy data from which relia- 
ble trends can be inferred for a wide range of model systems. 
The second and main objective is to achieve a detailed under- 
standing of the nature of halogen bonds: how they resemble, 
and also how they differ from, the better understood hydrogen 
bonds in terms of their electronic structure and bonding 
mechanism. 

To this end, we first explore how the geometries and ener- 
gies of our model complexes DX-A~ and, for comparison, 
DH-A~ vary as either the halogen- or hydrogen-bond donating 
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atom (D), or the halogen- or hydrogen-bond accepting atom 
(A) is varied from F to CI, Br and I. In this way, we arrive at 
a set of consistent data for a large range of halogen- and hy- 
drogen-bonded complexes. Next, to understand the origin of 
the computed trends, we carry out activation strain analyses 1251 
of the bond formation reaction. The activation strain model of 
chemical reactivity 1251 is a fragment-based approach to under- 
stand the energy profile of a chemical process in terms of the 
original reactants: the strain energy associated with their geo- 
metrical deformation and their mutual interaction along the re- 
action coordinate (see below). The interaction energy and the 
underlying bonding mechanism are furthermore analyzed in 
the context of quantitative Kohn-Sham molecular orbital (MO) 
theory in combination with an energy decomposition analysis 
(EDA). 126,271 Our explorations and analyses augment earlier pio- 
neering studies 128 " 311 through the large variety in our halogen- 
and hydrogen-bonded model complexes and the systematic 
and in-depth analyses along the entire reaction profile for each 
of the complexation reactions. 
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Theoretical Methods 

Computational details 

All calculations were carried out using the Amsterdam density 
functional (ADF) program developed by Baerends and co-work- 
ers. 132 " 341 The numerical integration was performed using the 
procedure developed by te Velde et al. [351 The MOs were ex- 
panded in a large uncontracted set of Slater-type orbitals 
(STOs) containing diffuse functions, TZ2P (no Gaussian func- 
tions are involved). The TZ2P basis set 1361 is of triple-^ quality 
for all atoms and has been augmented with two sets of polari- 
zation functions, that is, 2p and 3d on H, 3d and 4f on F and 
CI, 4d and 4f on Br, and 5d and 4f on I. The core shells of the 
halogen atoms (1s for F, 1s2s2p for CI, up to 3p for Br and up 
to 4p for I) were treated by the frozen-core approximation. An 
auxiliary set of s, p, d, f and g STOs was used to fit the molecu- 
lar density and to represent the Coulomb and exchange po- 
tentials accurately in each self-consistent field (SCF) cycle. 

Equilibrium structures were obtained by optimizations using 
analytical gradient techniques. 1371 Geometries and energies 
were calculated at the BP86 level of the generalized gradient 
approximation (GGA); exchange is described by Slater's Xa po- 
tential, 1381 with nonlocal corrections due to Becke 1391 added self- 
consistently, and correlation is treated in the Vosko-Wilk-Nusair 
(VWN) parameterization 1401 with nonlocal corrections due to 
Perdew 1411 added, again, self-consistently (BP86). [421 Scalar rela- 
tivistic effects were accounted for using the zeroth-order regu- 
lar approximation (ZORA). [431 Energy minima have been verified 
to be equilibrium structures through vibrational analysis. 1441 All 
minima were found to have zero imaginary frequencies. 

Throughout this paper, we focus on the electronic energies 
of the molecular systems. However, enthalpies at 298.15 K and 
1 atm are calculated using standard statistical mechanics rela- 
tionships as well. 1451 The thermodynamic effects were found to 
have only a small influence on the energies and do not alter 
the trends. For clarity, these results are, therefore, not dis- 
cussed but are included in the Supporting Information. 

Analysis of the bonding mechanism 

Insight into the bonding mechanism is obtained through acti- 
vation strain analyses of the various hydrogen- and halogen- 
bond formation reactions. These complexation reactions are 
computationally modeled by decreasing the distance between 
A" and the DH or DX fragment, and simultaneously increasing 
the D-H or D-X bond length. The DH-A~ or DX-A~ distance 
is decreased from an initial value of 1.8 times 1461 the equilibri- 
um bond length in the corresponding HA or XA molecule to 
the actual bond length value in the hydrogen- or halogen- 
bonded complex (r H ... A - or r x ... A -), while the DH or DX fragment 
is stretched from its equilibrium geometry to the geometry it 
acquires in the hydrogen- or halogen-bonded complex. Thus, 
each analysis starts from an optimized DH or DX molecule and 
a halide at a relatively large distance, which is then linearly 
transformed to the optimized hydrogen- or halogen-bonded 
complex. 



These complexation reactions are analyzed using the activa- 
tion strain model. The activation strain model of chemical reac- 
tivity 1251 is a fragment-based approach to understand the 
energy profile of a chemical process in terms of the original re- 
actants. Thus, the potential energy surface AE(t,) is decom- 
posed along the reaction coordinate t, (or just at one point 
along £,) into the strain energy AE strain (£,), which is associated 
with the geometrical deformation of the individual reactants 
as the process takes place, plus the actual interaction energy 
AE mt (Q between the deformed reactants [Eq. (1)]. 

A£(5)=A£ -nh (Q+A£ ht (Q (1) 

In the equilibrium geometry, that is, for t, = t, eq , this yields an 
expression for the bond energy A£(^ eq ) = A£= A£ strain + AE int . 
The PyFrag program was used to facilitate the analyses along 
the reaction coordinate t, of the bond formation processes. 1471 
The interaction energy AE int (Q between the deformed reac- 
tants is further analyzed in the conceptual framework provided 
by the quantitative Kohn-Sham MO model. 1261 To this end, it is 
decomposed in three physically meaningful terms [Eq. (2)] 
using a quantitative energy decomposition scheme developed 
by Ziegler and Rauk. 1271 

A£ int (0 = Al/ ebtat (g + AE Pauli (y + Af oi (0 (2) 

The term Al4i sta t corresponds to the classical Coulomb inter- 
action between the unperturbed charge distributions of the 
deformed reactants and is usually attractive. The Pauli repul- 
sion energy (A£ Pa J comprises the destabilizing interactions be- 
tween occupied orbitals of the reactants and is responsible for 
steric repulsion. The orbital-interaction energy (A£J accounts 
for charge transfer, that is, the interaction between occupied 
orbitals of one fragment with unoccupied orbitals of the other 
fragment, including the interactions of the highest occupied 
and lowest unoccupied MOs (HOMO-LUMO), and polarization, 
that is, empty-occupied orbital mixing on one fragment, due 
to the presence of another fragment. Since the Kohn-Sham 
MO method of DFT in principle yields exact energies, and 
rather accurate energies in practice, with the available density 
functionals for exchange and correlation, we have the special 
situation that an MO method, in principle, completely accounts 
for the bonding energy. 

The electron density distribution is analyzed using the Voro- 
noi deformation density (VDD) method for computing atomic 
charges. 1481 The VDD atomic charge on atom A (Q™) is com- 
puted as the (numerical) integral of the deformation density in 
the volume of the Voronoi cell of atom A [Eq. (3)]. The Voronoi 
cell of atom A is defined as the compartment of space bound- 
ed by the bond midplanes on and perpendicular to all bond 
axes between nucleus A and its neighboring nuclei. 

0f° = - J (p(r)-£ B A>(r))dr (3) 

Voronoi 
cell of A 

Here, p(r) is the electron density of the molecule and E^o^r) 
the superposition of atomic densities p B of a fictitious promole- 
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cule without chemical interactions where all atoms are consid- 
ered neutral. The interpretation of the VDD charge Q A DD is 
rather straightforward and transparent: instead of measuring 
the amount of charge associated with a particular atom A, 
Q A DD directly monitors how much charge flows out of 
(Q A DD >0) or into (Q A DD <0) the Voronoi cell of atom A due to 
chemical interactions. 



Results and Discussion 

Hydrogen-bond strength and structure 

The results of our ZORA-BP86/TZ2P calculations are shown in 
Table 1 for a representative selection of hydrogen-, fluorine- 
and iodine-bonded model complexes DX-A", covering D = F 
and I, X = H, F and I, as well as A" = F - , CI", Br and I" (full 
data on all model systems can be found in Tables SI -SI 0 in 
the Supporting Information). In the first place, we note that all 
model reactions are associated with single-well potential 
energy surfaces (PES), that is, there is no separate energy 
minima for DX-A - and D-XA - . In the case where D = A, D xh 
symmetric complexes with equal bond distances r D _ x = r x ... A - are 
formed (see Table 1). Furthermore, bond enthalpies at 298 K 
(AH 298 ) show the same trends as the electronic bond energies 



Table 1. Bond lengths [A] and energies relative to reactants [kcalmol 1 ] 
of the hydrogen-, fluorine- and iodine-bonded complexes. 1 * 

DX-A- r D _ x Ar D _ x Ibl r x _ k - Ar x ... A - [cl AE AE strain AE int BDE D _ x [d] 


FH-F" 


1.159 


0.226 


1.159 


0.226 


-53.0 


19.7 


-72.8 


143.5 


FH-cr 


1.012 


0.079 


1.843 


0.550 


-26.6 


3.3 


-29.8 


143.5 


FH-Br 


0.994 


0.061 


2.058 


0.625 


-21.9 


2.0 


-23.9 


143.5 


FH-r 


0.982 


0.049 


2.319 


0.694 


-18.1 


1.3 


-19.4 


143.5 


IH -F 


2.319 


0.694 


0.982 


0.049 


-80.6 


40.9 


-121.4 


81.7 


iH-cr 


2.191 


0.566 


1.423 


0.130 


-38.0 


31.7 


-69.6 


81.7 


IH-Br 


2.057 


0.432 


1.642 


0.209 


-28.6 


21.8 


-50.3 


81.7 


iH-r 


1.941 


0.316 


1.941 


0.316 


-21.8 


13.6 


-35.4 


81.7 


FF-F" 


1.755 


0.335 


1.755 


0.335 


-51.5 


23.5 


-75.0 


50.1 


FF-CI" 


1.864 


0.444 


1.965 


0.301 


-43.3 


34.2 


-77.5 


50.1 


FF-Br 


1.902 


0.482 


2.049 


0.253 


-44.0 


37.7 


-81.7 


50.1 


FF-r 


1.993 


0.573 


2.126 


0.181 


-48.4 


46.0 


-94.3 


50.1 


IF- -F 


2.126 


0.181 


1.993 


0.573 


-23.9 


5.9 


-29.8 


75.3 


iF-cr 


2.158 


0.213 


2.294 


0.630 


-14.5 


7.8 


-22.3 


75.3 


IF— Br 


2.200 


0.255 


2.335 


0.539 


-14.5 


10.5 


-25.0 


75.3 


iF-r 


2.324 


0.379 


2.324 


0.379 


-16.9 


19.3 


-36.2 


75.3 


Fl -F 


2.129 


0.184 


2.129 


0.184 


-75.0 


6.1 


-81.1 


75.3 


Fi-cr 


2.124 


0.179 


2.620 


0.268 


-49.8 


5.8 


-55.6 


75.3 


Fl— Br 


2.126 


0.181 


2.781 


0.275 


-45.1 


5.9 


-51.0 


75.3 


Fi-r 


2.132 


0.187 


2.977 


0.277 


-41.9 


6.3 


-48.2 


75.3 


1 1---R 


2.977 


0.277 


2.132 


0.187 


-69.0 


6.4 


-75.4 


49.0 


II— cr 


2.971 


0.271 


2.632 


0.280 


-44.2 


6.1 


-50.3 


49.0 


ll-Br 


2.976 


0.276 


2.795 


0.289 


-39.9 


6.4 


-46.3 


49.0 


II - 1 


2.991 


0.291 


2.991 


0.291 


-37.4 


6.9 


-44.3 


49.0 


[a] Computed at ZORA-BP86/TZ2P. [b] Stretch of the D-X fragment rela- 
tive to the optimized DX molecule, [c] Change in X— A distance compared 
to the situation in the optimized XA molecule, [d] BDE D _ X is the homolytic 
bond dissociation energy of the D— X bond without ZPE [kcalmol -1 ] 



A£ with differences in the order of one kcalmol -1 or less (see 
Tables S1-S5 in the Supporting Information). Thus, for clarity 
and a straightforward connection with our activation strain 
analyses, our discussion is based on the trends in electronic 
bond energies AE 

For the hydrogen-bonded DH-A - complexes, we find that 
as we vary the hydrogen-bond accepting halide A - from F - to 
T, the hydrogen-bond strength is weakened, r H ... A - becomes 
longer, and the D— H bond becomes less elongated from its 
equilibrium value in an isolated DH molecule (Ar D _ H = 
r D _ HA — r D _ H ). The opposite trend emerges as we vary the hydro- 
gen-bond donating atom D in DH down group 17 (F - to I"). 
Thus, along the hydrogen halides from FH to IH, the hydro- 
gen-bond strength is reinforced, the r H ... A - bond length be- 
comes smaller, and the bond stretch Ar D _ H increases. 

For example, from FH-F - to FH-I - , the hydrogen-bond 
strength (AE) is weakened from —53 to —18 kcalmol -1 , while 
the r H ... A - value increases from 1.159 to 2.319 A, and the value 
of the bond stretch Ar D _ H is reduced from 0.226 to 0.049 A 
(see Table 1). This trend correlates with a systematic weakening 
of the halide's proton affinity (PA) value of 373 kcal mol -1 for 
F - to 316 kcalmol -1 for I - (PA values at ZORA-BP86/QZ4P 
taken from Ref. [49]). The effect is even more pronounced in 
the series from IH-F" to IH-I" along which AE weakens from 
a value of —81 to —22 kcal mol -1 , the r H ... A - value increases 
from 0.982 to 1.941 A, and the Ar D _ H value is reduced from 
0.694 to 0.316 A. Note that, on the other hand, from FH-F - to 
IH-F - , A£ is strengthened from a value of —53 to —81 kcal 
mol -1 , while the value of r H ... A - decreases from 1.159 to 0.982 A, 
and the stretch value Ar c _ H is increased from 0.226 to 0.694 A. 
The higher extent of deformation in the more strongly hydro- 
gen-bonded complexes is also reflected by a more destabiliz- 
ing strain energy (A£ strain ; see Table 1). This trend correlates 
with a systematic weakening of the halogen-hydrogen bond 
from a homolytic bond dissociation energy (BDE) value of 
144 kcalmol -1 in FH to 82 kcalmol -1 in IH (see Table 2). Fur- 
thermore, note that the bond distance r H ... A - in DH-A - is in all 
cases longer than that of r H _ A in the diatomic HA molecule, as 
revealed by the corresponding difference in bond distances 
Ar H ... A - = r H ... A -— r H _ A . This difference Ar H ... A - increases in value 
from 0.226 A in FH-F - to 0.694 A in FH— I and from 0.049 A in 
IH-F - to 0.316 A in IH-I - . 

We conclude that the DH-A - hydrogen bond becomes 
stronger and relatively shorter, while the D— H bond becomes 
more elongated in the complex, as the A - anion is a stronger 
base and/or the D— H bond is weaker. 

Halogen-bond strength and structure 

In part, the halogen bonds display similar trends to the hydro- 
gen bonds, but there are also striking differences. In general 
and in agreement with ab initio results, the fluorine bonds are 
the weakest and the iodine bonds the strongest halogen 
bonds. 121 ' 23 ' 241 The heavier DX-A - halogen bonds (i.e., X = CI, Br 
and I) become weaker and longer as the accepting halide (A - ) 
varies from F - to I - , similar to the corresponding hydrogen 
bonds. In the case of the iodine-bonded complexes Dl-A", for 
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Table 2. Geometry, stability, and electronic structure of DH and DX mole- 
cules. 131 

D-X r D _ x BDE Q™ s(a) e(o*) sW e(at*) 



9.74 
8.04 
7.63 
7.03 

8.04 
7.37 
7.13 
6.78 

7.63 
7.13 
6.93 
6.62 

7.03 
6.78 
6.62 
6.39 

[a] Computed at ZORA-BP86/TZ2P; r D _ x = D-X distance [A]; BDE = homo- 
lytic bond dissociation energy without ZPE [kcalmol" 1 ]; Q™ = VDD 
charge on atom X [au]; t = orbital energy [eV]. 



example, A£ weakens from a value of around — JOkcalmoT 1 
for A~ = F~ to around -40 kcalmol" 1 for X" = l" (see Table 1). 
However, the fluorine bonds DF-A" display a more complex 
dependency of A£ upon variation of the accepting halide A". 
From A~ = F~ to CI", the fluorine bond strength still weakens, 
similar to the situation for the hydrogen bonds and the heavier 
halogen bonds. However, thereafter, along A~ = CI~, Br and 
r, the fluorine-bond strength no longer continues to weaken 
but instead becomes stronger. This is most clearly seen in the 
series constituted by the complex FF-A - between a fluorine 
molecule and a halide ion. Here, A£ for the fluorine-bond 
strength varies along A~ = F~, CI", Br and \" with values of 
-52, -43, -44 and -48 kcalmol" 1 , respectively (see Table 1). 

Interestingly, variation of the donating atom D has opposite 
effects on the halogen bonds DX-A" and hydrogen bonds 
DH-A". All halogen bonds studied here become weaker and 
longer as D is varied from F to I (see Table 1), whereas the hy- 
drogen bonds were found to become stronger and shorter 
along this series (see above). For example, along the series 
from FF-F" to IF-F", the fluorine-bond strength weakens from 
a AE value of —52 to only — 24kcalmor\ the fluorine-bond 
distance r x _ A - increases in value from 1.755 to 1.993 A, and the 
stretch Ar D _ x decreases in value from 0.335 to 0.181 A. 

Bond analyses with variation of A 

Our analyses show that the weakening of hydrogen bonds 
DH-A" and of heavier halogen bonds DX-A" (X = CI, Br, I), as 



the accepting group varies from A~ = F~ to I", is directly relat- 
ed to the concomitant reduction in electron-donating capacity 
of the np-type HOMO of the A" halide. The hydrogen and hal- 
ogen bonds appear to have an electrostatic component 
(Al/ e[stat ) and a covalent component (AfJ stemming mainly 
from the HOMO-LUMO interaction between the occupied 
halide np atomic orbital (AO) and the D— H or D— X antibonding 
a* acceptor orbital, shown schematically in Figure 1. Both 
bonding components, AV ekat and A£ oi , are weakened as the 
halide HOMO becomes more diffuse and effectively lower in 
energy 1501 from A~ = F~ to r (see Tables S6, S8-S10 in the Sup- 
porting Information). Consequently, the interaction energy 
(AE int ) and, thus, the net hydrogen- or halogen-bond strength 
A£ becomes less stabilizing along A~ = F~ to I" (see Table 1 
and Tables S1 and S3-S5 in the Supporting Information). 

The key to understanding why fluorine bonds DF-A" show 
a more complex, partially opposite trend (i.e., the expected 
weakening from A~ = F~ to CI" but thereafter a strengthening 
along A" = CI", Br" and l~) is contained in the counteracting ef- 
fects evolving from D— F bond stretching induced in the dia- 
tomic DF molecule as it interacts with the halide A". Interest- 
ingly, activation strain analyses reveal that from an early until 
a relatively advanced stage of the complexation reaction, for 
a given point along the reaction coordinate we indeed re- 
cover the original trend in interactions, namely, that AE int © 
weakens from A~ = F~ to I". This can be nicely seen in Figure 2 
which, for six representative series, shows the activation strain 
diagrams along the entire reaction coordinate t, projected 
onto the stretch Ar D _ x of the complexation reaction between 
a DX molecule approaching the halogen-bond accepting A" 
(see section, Theoretical Methods). Each of the six activation 
strain diagrams in Figure 2 refers to one particular DH or DX 
molecule forming a hydrogen or halogen bond with A" = F~, 
CI", Br" and l~. The strain-energy curves (A£ strain ) within each 
of these subgraphs coincide because they refer to the same di- 
atomic molecule being stretched as the complexation reaction 
progresses. Consequently, the trend A~ = F~ to I" in the total 
DH-A~ and DX-A~ energy profiles A£ in each subgraph is di- 
rectly determined by the trend in the corresponding interac- 
tion-energy curves (A£ int ). Also, as can be seen in Figure 2, the 
A£ in , curve appears to be most stabilizing for A~ = F~ and then 
weakens along CI", Br" and T, for any given diatomic mole- 
cule DH or DX, including all fluorine-bonded molecules DF. 

In other words, fluorine bonds DF-A" would also show 
a weakening in interaction A£ in , from A" = F" to \~, as the hy- 
drogen bonds and all other halogen bonds, if it were not for 
the increasingly stretched D— F bond in the fluorine-bond-do- 
nating diatomic molecule (see Table 1 and Figure 2). This struc- 
tural phenomenon is promoted by a combination of factors: 
1) a weak D— X bond that is easily stretched; 2) a strong inter- 
action with an approaching halide A"; and importantly, 3) a DX 
o* acceptor orbital that quickly drops in energy as the D— X 
bond elongates (see Figure 1). The latter generates a driving 
force for D— X stretching in DX-A" because it enhances the or- 
bital interactions and thus A£ in , (see Figures 1 and 2). Indeed, 
D— X stretching is most pronounced if this bond in the diatom- 
ic fragment is weaker, that is, for the weaker halogen-hydro- 



F-H 


0.933 


143.5 


0.20 


-13.57 


-0.72 


-9.78 


Cl-H 


1.293 


107.5 


0.10 


-11.79 


-0.97 


-8.05 


Br— H 


1.433 


94.6 


0.07 


— 11.18 


— 1.42 


—7.51 


l-H 


1.625 


81.7 


0.05 


-10.31 


-1.88 


-6.91 


F-F 


1.420 


50.1 


0.00 


-15.61 


-6.17 


-13.05 


Cl-F 


1.664 


69.2 


-0.07 


-13.61 


-4.86 


-11.66 


Br— F 


1 .796 


69.8 


— 0.1 1 


— 1 2.86 


—5.04 


— 1 1 .01 


l-F 


1.945 


75.3 


-0.13 


-11.95 


-4.86 


-10.49 


F-CI 


1.664 


69.2 


0.07 


-13.61 


-4.86 


-11.66 


CI-CI 


2.023 


62.0 


0.00 


-11.93 


-4.51 


-9.89 


Br-CI 


2.173 


58.8 


-0.03 


-11.38 


-4.71 


-9.36 


l-CI 


2.352 


57.9 


-0.08 


-10.72 


-4.67 


-8.93 


F-Br 


1.796 


69.8 


0.11 


-12.86 


-5.04 


-11.01 


Cl-Br 


2.173 


58.8 


0.03 


-11.38 


-4.71 


-9.36 


Br-Br 


2.321 


55.0 


0.00 


-10.88 


-4.82 


-8.86 


I Br 


2.506 


53.0 


-0.06 


-10.26 


-4.73 


-8.41 


F-l 


1.945 


75.3 


0.13 


-11.95 


-4.86 


-10.49 


Cl-I 


2.352 


57.9 


0.08 


-10.72 


-4.67 


-8.93 


Br-I 


2.506 


53.0 


0.06 


-10.26 


-4.73 


-8.41 


l-l 


2.700 


49.0 


0.00 


-9.68 


-4.65 


-7.92 



ChemistryOpen 2012, 1,96-105 © 2012 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim WWW.chemistryopen.org 99 




ChemistryOPEN 

L. P. Wolters, F. M. Bickelhaupt 

a) b) 

t>>Q -< no >C \ n 

n ^ HHHh P ^ 00 ^ , *^^fy 

Q • a -ff -ff 

DH DH—A" A" DX DX—A - A - 

c) _ d) 

icoc p ** / p ^ *+- p 8*~ 

a O ■ v.* , >-ff s Q 

D- D-H H" D - D-A A* 

Figure 1. Simplified orbital-interaction diagrams for a) hydrogen-bonded complexes DH— A~, b) halogen-bonded complexes DX-A", c) hydrogen halides D— H, 
and d) dihalogens D— X, as they emerge from our quantitative Kohn-Sham MO analyses. 



gen bonds (D— X = l— H) and the weaker halogen-halogen 
bonds (D-X=F-F; see Table 1). In the latter, it is able to affect 
the trend in overall bond strength A£. The D— F stretching in 
fluorine-bonded complexes is most pronounced in the FF-A - 
series, along which the F— F stretch Ar D _ x increases from 
a value of 0.3 via 0.4 and 0.5 to 0.6 A. This further stretch is 
able to induce the reversal of the trend in bond strength A£ 
along the equilibrium structures FF— CI - , FF-Br" and FF— I - (see 
Table 1). 

Thus, fluorine-bond analyses in the DF-A - equilibrium geo- 
metries show that in most cases the interaction energy (A£ int ) 
between the stretched D— F molecule and the halide A", as 
well as its components AV e | Stat and A£ oi , become more stabiliz- 
ing along the entire series A - = F - to I - , that is, already from 
F - to Cr (see Table 1, S2 and S7 in the Supporting Informa- 
tion). This is indeed most pronouncedly so in the series FF-A", 
due to the F— F bond in the DX fragment being relatively 
weak. Along the series FF-F", FF— CI - , FF-Br and FF— I - , A£ in , 
increases in strength from a value of —75 to —78, —82, and 
—94 kcal mor 1 , respectively. For comparison, along the corre- 
sponding series with the much stronger F— I bond in the DX 
fragment, that is, Fl—F , FI-CI - , Fl-Br" and Fl-I , the A£ int 
weakens from —81 to —56, —51, and —48 kcal mor 1 . 

The overall bond strength A£ along the fluorine-bonded 
series shows the aforementioned initial weakening followed by 
a strengthening, because the D— F stretching and the concomi- 
tant strain energy (A£ strain ) becomes more destabilizing along 



the series and, from A=F - to CI - , dominates the strengthen- 
ing in A£ int (see Table 1). 

We conclude that, in general, hydrogen bonds DH-A - and 
halogen bonds DX-A" become weaker along A - = F - to I" be- 
cause the larger radii and lower np AO energies of the halides 
lead to weaker electrostatic attraction and weaker orbital inter- 
actions. Interestingly, for the same reason, F is the halide with 
the strongest gas-phase basicity, the strongest alkyl cation af- 
finity and the lowest barrier for S N 2 reactions with halome- 
thanes. 149-511 The trend in DF-A - fluorine-bond strength is par- 
tially inverted, that is, A£ becomes more stabilizing along A" = 
CI", Br - and I - because of a more subtle interplay of factors. 
Notably, a significant stretching of the relatively weak D— F 
bond in the DF-A - equilibrium structures lowers the DF a* ac- 
ceptor orbital and thus amplifies the donor-acceptor orbital in- 
teractions, for example, along FF-CI - , FF-Br and FF---I . 

Bond analyses with variation of D 

We recall that for the hydrogen bonds DH-A - , a heavier do- 
nating halogen D results in a stronger bond, whereas the same 
variation in D weakens the halogen bonds DX-A - (see Table 1 
and S1-S5 in the Supporting Information). In both cases, the 
trend in bond strength A£ is determined by the interaction 
energy A£ int . For example, from FH-F - to IH-F - , A£ int is 
strengthened from a value of —71 to —121 kcalmol -1 , whereas 
from FI---F to II — R it is weakened from a value of —75 to 
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Figure 2. Activation strain analyses along the reaction coordinate for DX + A~ complexation as a function of A~ = F~, CT, Br and I - , projected onto the D— X 
stretch Ar DX for a) hydrogen bonds, b) fluorine bonds and c) iodine bonds, with donating groups D = F (left) and D = l (right). Energy profiles Af (solid lines) 
are decomposed into strain energy Af strain (dashed lines above A£=0) and interaction energy A£ int (dashed lines below A£=0). 



—29 kcalmor 1 (see Table 1). The strain energy (A£ strain ) is not 
negligible, but it does not alter the trend set by AE int . Our acti- 
vation strain analyses explain the above differences between 
hydrogen and halogen bonds, but they also confirm once 
more that both are very similar in nature (see Figure 3). 

Starting with some general observations, we find that for hy- 
drogen as well as halogen bonds, the strain-energy (AE strain ) 
curves are most unfavorable when D = F and gradually 
become less destabilizing as the donating atom is varied along 
D = F, CI, Br and I (see Figure 3). Furthermore, we find that for 
all DH-A~ and DX-A~ complexes, the interaction-energy 
(A£ int ) curves become less stabilizing along D = F, CI, Br and I. 
The resulting energy profiles and, therefore, the stability and 



geometric properties of the complexes DH-A~ and DX-A~ 
depend on the balance between the A£ strain and A£ int curves, 
which we first discuss individually. 

The slope and shape of the A£ strain curves is of course direct- 
ly related to the D— X bond strength of the diatomic fragment, 
which in general becomes stronger as the polarity across the 
D— H or D— X bond increases (see Table 2). This is a well-known 
and understood phenomenon. 1521 From FH to IH, the halogen- 
hydrogen bond strength decreases significantly from a value 
of 143 to 82 kcalmor 1 (Table 2). The corresponding halogen- 
halogen bonds are all much weaker, and variations in the ho- 
molytic BDE are also much smaller. From FF to IF, the bond 
strength increases from 50 kcalmor 1 to 75 kcalmor 1 , while for 
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Figure 3. Activation strain analyses along the reaction coordinate for DX + A~ complexation as function of D = F, CI, Br and I, projected onto the D— X stretch 
Ar DX for a) hydrogen bonds, b) fluorine bonds and c) iodine bonds, with accepting groups A~ = F~ (left) and A~ = T (right). Energy profiles AE (solid lines) are 
decomposed into strain energy Af strai „ (dashed lines above Af =0) and interaction energy AE int (dashed lines below A£=0). 



the fragments DX, where X is CI, Br or I, the bond strength 
generally decreases from a value of around 70kcalmol~ 1 for 
FX to around SOkcalmor 1 for IX. Thus, for the hydrogen- 
bonded complexes, the A£ strain curves show a pronounced re- 
duction in slope from FH to IH, which, in the corresponding 
hydrogen-bonded complexes FH-A" to IH--A", translates into 
an increasing stretch Ar D _ H of the diatomic fragment. As the 
stretch Ar D _ H becomes larger from equilibrium structures 
FH-A" to IH--A", the A£ int curves have been able to descend 
further, to lower, more stabilizing energies. The final result is, 
thus, an increasing stability of the DH-A" complexes when the 
donating atom D is varied from F to I. 



For the halogen bonds, the A£ strain curves are very similar 
and not decisive. The reason for the decreased stability of the 
DX-A" complexes upon the same variation of D from F to I is, 
therefore, that the A£ ln , curves descend more gradually to 
overall less stabilizing values. A£ in , becomes less stabilizing 
from FX-A" to IX-A" because of decreasing electrostatic at- 
tractions (AV e istat) an d, in some cases, also because of greater 
Pauli repulsions (Afp^; see Tables S6-S10 in the Supporting 
Information). Both of these effects are easily explained consid- 
ering the electronegativities of the halogens. Along the series 
FX to IX, the central atom X becomes relatively more electro- 
negative, which will lead to a greater negative charge on this 
central atom, thus, reducing the electrostatic attraction with 
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the anionic A", while concomitantly the occupied orbitals will 
have more X character, which in turn introduces stronger Pauli 
repulsion. 

Bond analyses with variation of X 

A more direct comparison of hydrogen and halogen bonds 
DX-A~ can be obtained by varying X along H, F, CI and I, while 
keeping the donating atom (D) and the accepting halide (A") 
constant (see Figure 4). We do this for four combinations of D 
and A", giving rise to four subgraphs in Figure 4. There ap- 
pears to be a regular trend of increasing strength from the flu- 
orine bonds to the iodine bonds. This trend derives again di- 
rectly from the electronegativity difference across the D— X 
bond of the diatomic fragment; from DF to Dl, the charge dis- 
tribution on the DX fragment is increasingly polarized towards 
D, away from A~ (see VDD atomic charges in Table 2), whereas 
the a* acceptor orbital achieves a higher amplitude on X (see 
Figure 1). This results in a strengthening of the halogen bond 
DX-A~ because of greater electrostatic attraction, less Pauli re- 
pulsion and more stabilizing donor-acceptor orbital interac- 
tions (see Tables S7-S10 in the Supporting Information). 

In analogy to the situation described above, hydrogen 
bonds might be expected to be much stronger than the halo- 
gen bonds due to the large and favorable polarization across 
the D— H bond leading to a partially positively charged hydro- 



gen atom in DH. For example, the VDD atomic charge on X in 
FH, FF and Fl amounts to +0.20, 0.00 and +0.13 au, respec- 
tively (see Table 2). The decomposition of the interaction 
energy into its components shows indeed a stronger electro- 
static attraction {AV elsat ) to the bonding energy in the case of 
the hydrogen bonds (compare Tables S6-S10 in the Support- 
ing Information). Note, however, that this does not imply that 
hydrogen bonds are always stronger than the corresponding 
halogen bonds, since in our model systems, the bonding 
mechanism is never purely, or even predominantly, electrostat- 
ic. The covalent or orbital-interaction term (A£ Q ) is relatively 
large and crucial for understanding the bonding in our model 
systems. For the hydrogen-bonded complexes DH-A+ the or- 
bital-interaction term accounts for 40 to 66% of the total 
bonding interactions {AV ektBt + AEJ. The stabilization due to 
this term results predominantly from a charge transfer from 
the np orbitals of the halide into the a* LUMO of the hydrogen 
halide (see Figure 1). For the halogen-bonded complexes 
DX-A", the contribution from the orbital-interaction term 
ranges from 43% for F1---F to as much as 97% for IF— F" at the 
other end of the spectrum. The larger covalent contribution, in 
the case of the halogen bonds, is the result of the low orbital 
energy of the empty dihalogen o* orbital (e.g., —0.7 eV for FH 
and —6.2 eV for FF; see Table 2), which directly translates into 
a stronger donor-acceptor orbital interaction with the halide 
np orbital (compare Tables S6-S10 in the Supporting Informa- 




Figure 4. Activation strain analyses along the reaction coordinate for DX + A~ complexation as function of X = H, F, CI and I, projected onto the D— X stretch 
Ar DX for bond donating groups a) D = F and b) D = l, and accepting groups A~ = F~ (left) and A~ = l~ (right). Energy profiles AE (solid lines) are decomposed 
into strain energy A£ strain (dashed lines above Af =0) and interaction energy AE im (dashed lines below A£=0). 
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tion). Note that, percentagewise, A£ oi in the halogen bonds ap- 
pears even larger because of the aforementioned, less favora- 
ble electrostatic attraction AV ektat . 

The nature of the strong hydrogen and halogen bonds dis- 
cussed in this work strongly resembles that of the weaker, neu- 
tral hydrogen and halogen bonds, although dispersion interac- 
tions become relatively more important in the latter. 12,19,24,311 
Preliminary results of dispersion-corrected ZORA-BP86-D3/TZ2P 
calculations on FI-FI (A£=-4.3 kcalmol 1 ), CICI-CICI (A£= 
-1.3 kcalmor 1 ) and 11-11 (A£= -6.6 kcal mor 1 ) show that the 
covalent component A£ oi amounts to 43-59%, whereas disper- 
sion contributes 2-17% to the total of all bonding interactions 
(A£ oi + AI/ e | Stat +dispersion interaction; see Theoretical Meth- 
ods and ref. [53] for the D3-dispersion correction as proposed 
by Grimme). The covalent contribution in these neutral model 
complexes stems from a donor-acceptor orbital interaction 
from an occupied n* orbital on one dihalogen fragment into 
the a* orbital of the other dihalogen fragment. 

We conclude that halogen bonds DH-A~ and hydrogen 
bonds DX-A~ have a very similar bonding mechanism consist- 
ing of both electrostatic and covalent contributions. The elec- 
trostatic attraction is less favorable in the halogen bonds due 
to a smaller and in some cases less favorably oriented polariza- 
tion across the dihalogen molecule DX. Nevertheless, halogen 
bonds can become stronger than hydrogen bonds because of 
a more stabilizing covalent component in the former. The 
reason is the lower orbital energy of the empty a* orbitals in 
dihalogen molecules DX leading to a stronger, more favorable 
donor-acceptor orbital interaction with the halide A~ np orbi- 
tal (see Table 2). 

Conclusion 

Halogen bonds in DX-A~ are very similar in nature to hydro- 
gen bonds in DH-A" (D, X, A = F, CI, Br, I): both have a sizeable 
covalent component stemming from HOMO-LUMO interac- 
tions between the np-type lone pair on the halogen- or hydro- 
gen-bond accepting fragment A~ and the D— X or D— H anti- 
bonding a* LUMO on the halogen- or hydrogen-bond donat- 
ing fragment DX or DH, respectively. Neither halogen bonds 
nor hydrogen bonds are, therefore, predominantly, let alone 
purely electrostatic phenomena. This follows from our bonding 
analyses based on relativistic density functional theory (DFT). 

Two characteristic differences between the halogen bonds 
DX-A~ and hydrogen bonds DH-A~ is that the former are 
generally associated with a weaker electrostatic attraction (di- 
halogens DX are less polar than hydrogen halides DH) and 
a significantly more stabilizing HOMO-LUMO interaction. The 
stronger orbital interaction derives from the lower energy of 
the halogen-halogen a* LUMO as compared with that of the 
much stronger halogen-hydrogen bond. Halogen bonds can 
be stronger but also weaker than the corresponding hydrogen 
bonds. 

Finally, hydrogen bonds DH-A~ and halogen bonds DX-A~ 
become weaker along A~ = F~ to I", because the electron-do- 
nating capability (and basicity, alkyl cation affinity, nucleophilic- 
ity) [49,511 of the halide decreases in this order. The trend in 



DF-A" fluorine-bond strength is partially inverted, that is, A£ 
becomes more stabilizing along A~ = CI ", Br and r, because 
of a more subtle interplay of factors, in which a significant 
stretching of the relatively weak D— F bond lowers the DF o* 
acceptor orbital and, thus, amplifies the donor-acceptor orbital 
interactions. 
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